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SUM-OF-EXPONENTIAL STATIONARY DENSITIES 
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Abstract. We give necessary and sufficient conditions for the stationary density of semi- 
martingale reflected Brownian motion in a wedge to be written as a finite sum of terms of 
exponential product form. Relying on geometric ideas reminiscent of the reflection princi- 
QO ' pie, we give an explicit formula for the density in such cases. 
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It is well-known that Brownian motion on the positive half-line with negative drift and 
reflection at zero has a stationary density which is exponential. One derivation of this 
fact uses time reversal to relate the distribution of the corresponding transitory process at 
time t in the driftless case to the distribution of the maximum of a standard Brownian 
motion over [0, t] , which can be found using the reflection principle. A Cameron-Martin- 
Girsanov change of measure then introduces the drift, and letting t — > oo gives the required 
' stationary distribution. In the (i-dimensional setting, analogous arguments show that in 

certain situations, semimartingale reflected Brownian motion (SRBM) in a polyhedral cone 
has a stationary density which can be written as a finite sum of terms of exponential product 
form — that is, terms of the form x i— > ae~^ x,x ^ for some o € 1, A £ We call such a density 
a sum of exponentials. The aim of this paper is to give necessary and sufficient conditions 
under which the stationary density of SRBM in a two-dimensional wedge can be written as 
a sum of exponentials. 

In order to motivate our study, we describe the above argument in a multidimensional 
setting in Section [2] below. The invariant measure for SRBM with special pushing directions 
at the boundary of a polyhedral cone then relates to exit probabilities for the corresponding 
free Brownian motion. For some cones, in analogy with the one-dimensional case, these 
exit probabilities can be obtained explicitly using the reflection principle. If the reflected 
process is positive recurrent, then upon differentiating we obtain a stationary density which 
can be written as a sum of exponentials. Since the latter can thus be viewed as a natural 
multidimensional extension of the exponential stationary density of one-dimensional SRBM, 
this raises the question under what circumstances the stationary distribution of a given 
SRBM in a polyhedral cone has a sum-of-exponential density. This paper shows that, for 
SRBM in a two-dimensional polyhedral cone, sum-of-exponential stationary densities arise 
well beyond those cases to which the above arguments can be applied. 

The stationary distribution of multidimensional SRBM in a polyhedral cone is only known 
in relatively few cases. It has a density consisting of a single exponential term under a skew- 
symmetry condition on the pushing directions on the faces [17\ I25j . We also mention an 
apparently isolated example due to Harrison [13] and work of Foschini [11] ; we discuss the 
latter in more detail below. In the absence of further explicit results, numerical techniques 
have been developed [6] and logarithmic tail asymptotics have been investigated [31 [9l Q3] . 

Before describing our results in more detail, we introduce some notation which is sum- 
marised in Figure [TJ Define the wedge as 

S = {x £ M 2 : < arg(x) < £}. 

Throughout this paper we write wq = (cos 9, sin 9) for Set. Given some < S, e,£ < vr, 
we set V\ = \\vi\\ws, i>2 = H^H^f-e- Let F\ and F2 be the two faces of the wedge. As is 
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Figure 1, The wedge S. 



customary, we normalise v\ and V2 such that (vi,m) = 1 and (^2,^2) = 1> where n\ and 
ri2 are the unit normal vectors on F\ and F%, respectively. Let some vector \i £ M 2 also be 
given, and write 6^ = arg(^u) E (— 7r, 7r]. A key role in this paper is played by a, which is 
introduced by Varadhan and Williams [22] as 




Our results concern the case a < 1, and it is well-known (e.g., [ 22 j 24j) that under this 
condition there exists a continuous semimartingale Markov process with properies that can 
intuitively be summarised as follows (we do not give the mathematically rigorous definition 
here, as it can be found in, e.g., [31 1101 123] ). The process behaves like a standard Brownian 
motion with drift — /z in the interior of the wedge S, at the boundary of the wedge it is 
pushed in some specified direction, and the time it spends at the vertex of the wedge has 
Lebesgue measure zero. The pushing directions are constant along each of the faces, and are 
given by V{ for the i-th face. We call this process SRBM in a wedge and seek its stationary 
distribution, which is absolutely continuous with respect to Lebesgue measure [5l 116]. 

We prove that the stationary density of this process is a sum of exponentials if and only 
if a = — I for some nonnegative integer t. Moreover, it follows from our results that the 
number of exponential terms is 2^ + 1 when a = —I. Note that this condition reduces to 
5 = it — e if £ = 0, i.e., the skew-symmetry condition for product forms |17j . In fact, we 
show that the stationary density can be written as an (I + 1) x {I + 1) determinant with a 
special structure. A corollary is that the density behaves near the origin as in the driftless 
case studied earlier by Williams [23] . 

We next discuss some work related to our additive generalisations of the product form 
exponential stationary density. To our knowledge, Foschini's study of symmetric SRBM 
in the wedge [11] was the first to consider sum-of-exponential stationary densities. Our 
paper is a continuation and extension of Foschini's work, in the sense that we couple key 
ideas from his algorithm with geometric arguments suggested by the reflection principle, and 
consequently obtain an explicit formula for the stationary density. A discrete-state space 
version of Foschini's algorithm, the compensation method, has been successfully applied 
to a variety of queueing problems [U [2]. Another body of work loosely related to the 
present paper studies (driftless) two-dimensional SRBM with a = — 1 [HI [19], motivated by 
a connection with Schramm-Loewner evolutions. 

This paper is organised as follows. Section [2] gives some background and discusses the 
special sum-of-exponential stationary densities mentioned above. Section contains our 
main result and gives a geometric construction of the stationary density. In Sections 2HH1 
we prove the main result. 
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2. Survival probabilities, reflection groups, and Weyl chambers 

This section describes a special class of sum-of-exponential stationary densities which can 
be obtained directly using time reversal. These densities motivate our main result, since 
they typify the general structure of a sum-of-exponential stationary density. 

Consider an SRBM as in Figure [H with drift — \i and 5 = e = £. Suppose that /i G S° 
(the interior of S) to ensure positive recurrence, and write II for the stationary measure of 
the SRBM. Time reversal gives 

IT({y G S : (y,ni) < (x,m), (y,n 2 ) < (x,n 2 )}) = P- X {T = oo), (1) 

where P y is the law of the corresponding free Brownian motion B starting in y and T is the 
first exit time from —S, i.e., T = inf{i > : B(t) —S}. To see this, we use the wedge — S 
to define a partial order < on R 2 (see for example |20| ) and then if B s = Bt- S — Bt, 

sup{B s : < s < t} - B t < x if and only if (2) 
sup{-x + B s :0<s<t} < 0, (3) 

where sup denotes the supremum with respect to the partial order. On the left hand side 
of © is the supremum over [0, t] of the free Brownian motion starting in —x, while on the 
left hand side of ((2|) we have an SRBM in S as defined above (starting in 0; for details see 
[20]). Applying Wiener measure to ([21) — (131) and letting t — > oo gives (P). 

Biane et al. [1] have recently shown that if the wedge angle is of the form £ = ir/m for 
some integer m > 2 then 



P_«(T = oo) = sgn (w)e-^ I - w ^ (4) 
weG 



for any x G S° and any /i G S°, where the sum is taken over the finite reflection group G 
with associated signature function sgn(-), as detailed below. We remark that this formula 
gives the probability that a standard Brownian motion with drift never exits a so-called 
Weyl chamber — in particular, it is not restricted to a two-dimensional setting. 

We next describe the index of summation in ([3]), i.e., the elements of the reflection group 
G. Throughout, we represent any element of M. 2 as a column vector. Write Rq for the 
reflection matrix across the line with argument 8, i.e., 



Ra 



cos 29 sin 29 
sin 29 - cos 29 



and pe for the rotation matrix over an angle 9, i.e., 



pe 



cos 9 — sin 9 
sin 9 cos 9 



The group G consists of the 2m matrices /, R%, p2£> p2£-^£> • • • 3 P2(m-i)£> P2(m-i)£-ft£- The 
signature of the matrices p2k£R£ is —1, while the signature of the others is +1. 

On combining ([T]) and (J3|) , we find that the stationary density of SRBM with £ = 5 = e = 
ir/m is proportional to 

sgn (w)(p, (I - w)v 2 ){fi, (I - w^e-M-vM . (5) 

weG 

This is a sum of 2m— 3 nonzero terms, since the prefactor vanishes when w = I , w = R%, and 
w = /02(m-i)f-% • I* i s the aim of this paper to put the explicit formula ((5J) into the context 
of more general SRBMs. Specifically, we show that every sum-of-exponential stationary 
density has a representation reminiscent of ©. 
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3. Main results 



This section presents our main result on sum-of-exponential stationary densities, and 
explains its geometric interpretation. Our result shows that sum-of-exponential stationary 
densities can be written as determinants with a special structure. The proofs are deferred 
to Sections SH3 

The matrices 

Rotfc = p25+2ki and Ref^ = P2<5+2(ifc-l)£-% (6) 
for k > play an important role in our main theorem. Write 

e £ = {6 G (f - e, S) : sm(6 - 26 - k£) + for k = 0, . . . , 21}. 

We know from Hobson and Rogers [18] and from Dupuis and Williams |10j that there exists 
a unique stationary distribution of the SRBM if £ — e < 9^ < 8. Let p^ be its density with 
respect to Lebesgue measure, and write e± for the vector (1,0). For integers j, we also define 

R through 



the function ttj : S 



irAx) 



{H, (I - Rot^v^e-^'V-* ^ - (/x, (J - Rei^vije-^'V-™^ 
(/i, (Refj — Rot,)ui) 



Observe that both p^ and tt^ depend on 8, e and £, which is suppressed in the notation. 

Theorem 1. The following are equivalent: 

(i) a = — I for some integer £ > 0; 

(ii) a = —£ for some integer I > and for all fj, with 6^ £ @£, the functions p^ and 7r^ 
are equal up to a multiplicative constant, where for x € S, 



TT r r . 



(X) 

(^,Rot ei) f 1 

(/i,Rot ei) 2 
(//,Rot ei) 
1 



7lf (a?) 

^Rotiei)^ 1 

(/j,Rotiei) 2 
(/j,,Rotiei) 
1 



7T^(x) 

(/x,Rot^ei) 



(^,Rot £ ei) 2 
(/i,Roteei) 
1 



(7) 



(iii) a < 1 and for some \x with £ — e < # M < 6, there exist K < oo, coefficients ai, . . . , ax, 
and vectors d\,..., dx such that p^ admits the representation 

A 



p^{x) 



i=i 



In part (ii) of this theorem, the condition 9^ £ prevents a certain degeneracy which 
is fundamental to the problem. Indeed, it guarantees the linear independence of the expo- 
nential functions in 7r^. We remark that the determinant in (ii) has a different form than 
the ones recently studied in connection with transition probabilities for certain Markov 
chains [7J. 

A straightforward calculation using Vandermonde matrices shows that the function ir^ 
defined by ([7]) may be written as 



y~] c k [(/i, (I - Rot fc )ui)e" 



(/x, (J - Ref fc )«i)e' 



-<jjL,(I-Be£ k )x) 



k=0 



where 



Ck 



-If 



k Ho<i<j<£;i,jytk(^ ( Rot i ~ Rotj)ei) 



(8) 



(9) 



(/i, (Reffc - Rotjs,)^) 

This representation for 7r^ is used later in the proof of Theorem[TJ note that each exponential 
term is characterised by a rotation or a reflection matrix, Roto, Ref i, Roti, . . . , Ref^, Rot^. 

Although there is not necessarily an underlying reflection group in the general setting of 
Theorem [TJ the rotation and reflection matrices in Theorem [T] suggest a connection with 



REFLECTED BROWNIAN MOTION IN A WEDGE 



5 




the reflection- group framework. We illustrate this for a = —2 (hence I = 2) in the leftmost 
diagram of Figure [21 where we depict the five points Rotox, Refix, Rotix, Ref2X, Rot2X for 
an arbitrarily chosen x E S. By construction, Rotox lies in the wedge p 2 sS, which we call the 
'initial wedge'. The other points Refox,Rotix, . . . are constructed by successive reflections, 
reminiscent of the orbit of Rotox under the action of a reflection group. In particular 

• Each point lies in one of the wedges constructed by rotating the initial wedge anti- 
clockwise over multiples of the wedge angle £. 

• Two points lying in adjacent wedges — that is, wedges which share a common bound- 
ary line — are reflections of each other in that line. 

Note that Refox, although not contributing to 7r^ since (/ — Refo)f i = 0, is also obtained 
by reflection from Rotox. We indicate this by a dashed line in Figure [2j 

If a = — £ then the last point lies in the wedge P2S+2£^S, which is the same as p~2eS. On 
comparing this with the initial wedge P25S, it transpires that the last wedge in the anti- 
clockwise construction given above is the first wedge in the following clockwise construction. 
For k > we introduce the matrices 

Rotfc = P-2k£-2e and Reffc = p_2(fc-l)£-2e-Ro- 

The rightmost diagram of Figure [2] illustrates a clockwise construction for the stationary 
density, starting with the wedge p~2eS and labelling the points using the matrices Rotfc and 
Reffc. The problem is exactly the same as in the leftmost diagram, so by uniqueness the 
corresponding sum-of-exponential densities must agree on S. We use this observation in the 
proof of Theorem [TJ 

We close this section by stating some properties of the function 7r M defined in ([7]), which 
play an important role in our proof of Theorem [TJ These properties are proved in Section [H 
We start with a result for the coefficients {c k } defined in ([9]). 

Lemma 1. Let —a = £ 6 {1, 2, . . .} and 6^ G &£. The coefficients Ck defined by ([9]) satisfy, 
for 1 < k < £, 

Ck{p, (I ~ Ref fc )wi)(|U, (I - Rot A ._ 1 )t) 2 } = c k -i(fi, {I - Ref k )v 2 )(p, (I - Rotfc_i)«i). (10) 

Using this lemma, it is readily checked that (0) is recovered from Theorem [JJ by setting 
5 = e = £ = 7r/m. Indeed, we then have {p, (I — Rotfc_i)ui) = {p, (I — Reffc)t?i), so that 
(fT0[) reduces to 

Chip, (I ~ Rotfc^)^} = c k _i{p, (I - Reffc)« 2 ). 
In this special case we may therefore set 

Ck = {p, (I ~ Rot k )v 2 ) = (p, {I - Ref k )v 2 ), 

and we obtain ©. 
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The limiting behaviour of 7r^ given in the following proposition should be compared with 
the invariant measure found by Williams [23] for reflected Brownian motion without drift. 
We abbreviate lim r ^o f(r)/g(r) = 1 by f(r) ~ g(r) as r — > 0. 

Proposition 1. If —a = £ G {1, 2, . . .}, then for any 9 G [0, £] and any // wni/i #^ G 0£, we 
aaue ^(rwg) ~ CV^ sin^+o') as r — > 0, where is a finite nonzero constant independent 
of r and 6. 

Our next result is that tt 11 defined by (J7|) does not change sign on S. Note that this 
resolves Conjecture 1 in Dai and Harrison [6] for the special class of SRBMs studied in this 
paper. 

Proposition 2. The function ir^ does not change sign on S. 

4. Properties of ir^ 

In this section, we prove the properties of tt^ claimed in Section El The proof of the main 
result, Theorem [H is deferred to Sections [5] and El 

4.1. Proof of Lemma [TJ We first divide (|10p by s'm(6^ — 6 — k£) sin(5 + (k — 1)£), which 
is nonzero as a consequence of the assumption on \x in conjunction with the identity 5 + e = 
7r — ££. Again using this identity, we find after some elementary trigonometry that (|10p is 
equivalent to, with u>t = Op — 25 — 

c fc sin(&;£) sin(t^ +fc ) = -c fc _i sin((£ + 1 - fc)£) sin(o; i; _ 1 ). 

To show that this holds for the defined in ([5]), we observe that 

(fi, (Reffc - Rx>tfc)«i) sin(o;2(fe-i)) = (m, (Ref fc _i - Rotfc_i)ui) sin(w 2 fc) 

and that 

(//, (Rotj - Rotj)ei) = -2sin((j - sin(o; i+ j). 
After some algebra we also find that 

sin(fc£)sin(^ +fc )sin(w 2 (fc-i)) ] [ (Rot, - Rot,-)ei) 

0<i<j<£;i,j^k 

= sin((£ - k + 1)£) sin(w fe „i) sm(u 2 k) ] | (m, (Rot; - Rotj)ei), 

Q<i<j<£;i,j^k-1 

and the claim follows. 

4.2. Proof of Proposition [TJ For simplicity we suppose that \\fj,\\ = 1. We first investigate 
the behaviour near zero of tt^(x), for which we rewrite e^' x ^ ir^ (x) using the determinantal 
representation ([7]). A key ingredient is the identity e rcos ^ = /o( r ) + 2 Y2n=i cos ( nr ])In( r ), 
where I n is the modified Bessel function of the first kind. Using this identity, after absorbing 
e^' x i into the first row, we rewrite the elements on this row as 

\\x\\(ii,Kotjw e ) _ p \\x\\(^KeijW B ) 

<„,, ?w = J*>W- W - M|W ( ,, (Itot; - M>l) <") 



e 



I (||x||) + 2 J2 ^n(cos(u; 2j - 0))I„(||z||) 

n=l 

oo 

/o(|M|) + 2 £ T n (cos(u; 2j - 0))I„(||x||) 

n=l 

oo 

- 2(fi, (I - Reij)vi) sin(n6')C/ n _i(cos(L<j 2 j))/„ 



n=l 
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where we again set uik = 6^ — 25 — k£, and T n and U n are the Chebyshev polynomials of the 
first and second kind, respectively. 

In conjunction with some trigonometry, the above reasoning shows that 



>^<(x)=/ (IM|) + 



sin(<5) 



J>i,„(0)I„(||x||), 



n=l 



where hj in (9) is defined as 

- sin(n9+S)U n (cos(u}2j))-(fi,vi/\\vx\\) sm(nO)U n ^ 1 (cos{u 2 j)) + - sin(n^-5)C/ n _ 2 (cos(a;2j))- 

(We use the convention J7_i(x) = 0.) Therefore, e^' x ^Tr fl (x) can be expanded in terms of 
modified Bessel functions of the first kind, and for n > 1 the coefficient in front of J n (||x||) 
is proportional to 



ho,»(0) h 2 , n (0) 

L^_i(cos(o;o)) C^_i(cos(w 2 )) c^_i(cos(w 4 )) 
C/ (cos(w )) C/ (cos(a;2)) £/o(cos(w 4 )) 



M ) 

LVi(cos(u^)) 
C/ (cos(o;2^)) 



(12) 



To see how this follows from (|7|), note that we may apply elementary determinantal oper- 
ations to replace a row with elements cos(i02j) m by L^cos^j)). The term /o(IMI) is not 
present in the expansion in view of the last row in (|12p with ones. 

The condition 6> M E 6^ guarantees that none of the cos(w2j) are equal, and we conclude 
that the coefficient of I n (||x||) vanishes for n < £ and that it is proportional to s'm(£0 + 5) 
for n = I. Since Iiif) ~ Cr^ for some constant C ^ as r — > 0, this yields ^^(rwg) ~ 
CV S in(^ + ( j). 

4.3. Proof of Proposition [2j We rely on the following auxiliary lemma, whose proof is 
inspired by some elementary symmetric-function theory. Alternatively, as communicated 
to us by Sean Meyn, Proposition [2] can be proved using a continuous-space analogue of 
Theorem 1 of Foster [12]; this can be derived using the general theory of Markov processes. 



Lemma 2. Let £ G {1,2, . . .} and £ G 



For any y > 0, the sign of the determinant 



e Coy 


e tiy . 


■ e^ y 




ct 1 ■ 


■ c l ~ x 


Co 


Ci • 


■ Ce 


I— 1 


1 


1 



equals the sign of rio<i<j< JO ~ Cj]- 

Proof. The statement is a continuous analogue of the claim that S( n ,o,...,0) (C) i s nonneg- 
ative for £ > 0, where s\ is a symmetric polynomial known as a Schur polynomial (or, in 
this special complete homogeneous symmetric polynomial). 

By induction on £ one can show that the given determinant equals 



IT & 

o<i<j<e 







3 Co(«o-*-i) . . . e Ct( z i- z t-i)d 



dzp 



0=z- 1 <z <...<z e _ 1 <z e =y 



and the claim follows. 



□ 



We now prove Proposition [2J By the Maximum Principle (see Theorem 2.5 of [2T] for a 
suitable form), neither the minimum nor the maximum of tt^ over S is attained in the open 
set S°. We therefore investigate the boundary. 
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We first prove that tt^(x) — > as — ► oo by showing that, for any x G S, 

(H, (I - Rot k )x) > 0, k = 0, . . . , £, 
(/j,, (I - Ref k )x) > 0, k = !,...,£. 

Set 9 = avgx. For the claim involving Rotfc, we observe that 

(jji, (I - Rot k )x) = -2\\x\\ sin(<5 + fcf ) sin(^ - 9 - S - fcf ). 

Since 0<5 + k£<5 + ££ = ir — e < ir and we have 9 < £ and the stability condition 
£ — e < 9^ < 5, we obtain — tt < 9^ — 9 — S — k£ < 0. The same argument works for the 
claim involving Reffc, now relying on 

(jjt, (I - Rei k )x) = -2\\x\\ sm(5 + k£ - 9) sin(^ — 5 — k£) 

and the assumption k > 1. 

We next prove that the signs of ^(rwo) and ir^frw^) are equal and independent of r > 0. 
The equality of the signs follows from Proposition [T] after showing that they do not depend 
on r. From Lemma[2]with Q = (/i, Rot,ex) and we conclude this for ^{twq). Applying 
the same argument 'clockwise' shows that this also holds for ir^(rw^). 



5. The BAR and a PDE with boundary conditions 

This section prepares for the proof of Theorem [T] by relating the stationary density to a 
partial differential equation (PDE) with boundary conditions involving the pushing direc- 
tions. 



5.1. The BAR. Our proof of Theorem [T] requires the Basic Adjoint Relationship (BAR) 
as presented in following proposition, which is implied by Propositions 3 and 4 in [6]; see 
[5] for proofs. 

Proposition 3. Suppose that a < 1 and assume the existence and uniqueness of a stationary 
distribution for the SRBM. 

A nonzero finite measure v$ on S is proportional to the stationary distribution if and only 
if there exist finite measures v\ on F\ and V2 on Fi such that for any f G C?(S) 

[ Ia/-( M ,V/) duo+ f (v 1 ,Vf)du 1 + f {v 2 ,Vf}dv 2 = (BAR). 

JS l z J J Fx JF 2 

Let a and Oi be the Lebesgue measures on S and Fi, respectively. Write v* = 2n« — vi. 

Proposition 4. Let p £ C 2 (S) be nonnegative and integrable over S. 
If (BAR) is satisfied with dv$ = pda and dvi = p/2da{, then 

Ap + 2(n,Vp) = onS°, (13) 
(vl,Vp) + 2(fj,,m)p = onFf, (14) 
(i>2,Vp) + 2{fi,n 2 )p = onF 2 °. (15) 

Conversely, if [T3\) - [15\) hold and moreover p(0) = 0, then (BAR) is satisfied with dvo = 
pda and dvi = pj2doi. 

Proof. We may repeat the arguments in the proof of Lemma 7.1 of Harrison and 
Williams [17]. The additional assumption p(0) = ensures that (7.8) in [T7] automati- 
cally holds. □ 

The above proposition motivates investigating sum-of-exponential solutions to (|13|) - (|15|) . 
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5.2. Sum-of-exponential solutions to the PDE. In this subsection we study some prop- 
erties of sum-of-exponential solutions to the PDE (|13f) plus either boundary condition (I14p 
or (fT5l) . We use the following observation, due to Foschini [TU Sec. III. A]. 

Lemma 3. Let p be given by p(x) = £^Li cbie~( Ci,x ) for some k < oo ; aj 7^ 0, and q 7^ 
suc/i i/iaf q 7^ Cj «/ i 7^ j . 

i/p satisfies $13\) and ( fi^| ), f/ien /or eac/i z = 1, . . . , k precisely one of the following holds: 

(1) x 1 — > e~( Ci ' x } satisfies j!4\ ), or 

(2) there exists a unique j ^ i such that x 1— ► a^e - ^ 1 ' 2 ^ + aje~^ c i ,x > satisfies [Lfy and we 
have (ci,w ) = (cj,w ). 

By symmetry, Lemma [3] also holds when (I14D is replaced by (]15p . provided the condition 
(ci,u>o) = (cj,Wo) is replaced by (ci,w^) = {cj,w^}. The next two lemmas investigate the 
two scenarios of Lemma [3] in more detail; Lemma H) may be regarded as a generalisation of 
Theorem 6.1 of [T7] (modulo the discussion of the BAR in Section [6]). 

Lemma 4. Let p(x) = e~^ c,x > for some c 7^ 0. 

(1) Ifp satisfies JTJJ and {7J), then either p(x) = e -<^( J -«M)*> orp(x) = e -0*.(^-flo)*> . 

(2) Ifp satisfies and (Q3]) ; t/ien either p(x) = e -<M,U-p-2 e )^) orp ( x ) = e -<Ai,(/-fle)*> _ 

Proof. We only prove the first claim, the second being the clockwise analogue. The 
condition that p satisfies (|13p translates to \\c — p\\ = We may therefore write c = 

p — P-2-yP for some 7 = 7(/f). Next we substitute this in ([14")) . yielding —{pi, (I — P2i)v*) + 
2(n,ni) = 0, which we may rewrite as (p,, v% + P2^v\) = since v* = 2ni — ui. Using 
v\ = —RsP28V\, we get (p, (I — P2(-i-5))v\) = 0. This can only be the case if 7 = 6 mod ir 
or 7 = 9^ mod 7r. □ 

The next result investigates the second scenario of Lemma El 

Lemma 5. Let p be given by p(x) = a\e~^ c ' x ^ + o^e - ^' 3 ^ for some 01, a2 7^ 0, c, d 7^ 0, and 
c </. 

(1) // (c, u;o) = (d, uio) and p satisfies and |7^P , £/ien £/iere exists some 7 = "f(p) £ 
(0, 7r) such that p is proportional to p 7 defined by 

p 7 (x) = (p, (I - P2 1+ 2&)vi)e-^ 1 -^^ - (p, (I - P2 1+ 2sRo)vi)e~^ I -^ R ^ x l 

(2) // (c,w^) = {d,w^} and p satisfies ( TO)) and 173)) . £/ien i/iere exzsis some 7 = j(p) £ 
(0, 7r) such that p is proportional to p^ defined by 

Py(x) = (p, (I - p^^Me-^ 1 -"-^-^ - (p, (I - p^^eR^e-^ 1 -^-*^. 

Proof. Again we only prove the first claim. By linear independence both e~^ c,x < and 
e —(d,x) mus ^ sa tisfy (|13h individually. As in the proof of Lemma H) we may therefore write 
c = p — p-26-2jP for some 7 = 'y(p) E [0, 7r), so that (c, x) = (p,(I — P2y+2s)x). From 
(c, ^o) = (d, wo) and c^dwe conclude that (d, x) = (p, (I — p2^+2sRo)x). 

It remains to study a\ and 02, for which we use (fl4"l) . Since (c,wq) = (d, u>o) we obtain 
that on 

(«J,Vp(s)) = -[ai(/i, (/ - P2 1+ 2s)vl) + a 2 (p, {I - p2~ f +2sRo)v* 1 )]e- {c > x) . 

With v\ = 2n\ — v\ , we conclude that (fl4")) implies 

ai(/i, Ul + /0 27 +25Wl) + a2(M, «i + P2 1 +2sRqv{) = 0. 

The result follows after using v\ = —RqV\ and noting that (p, (I — p2- l +2sRo)vi) cannot 
be zero. □ 

We remark that Lemmas EHSl show that some structure from the particular example of 
Section [5] holds in general. Specifically, each exponent in a sum-of-exponential solution 
equals —{p, (I — M)x) for some reflection or rotation matrix M. 
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6. Proof of Theorem Q] 

To prove our main result, it suffices to show that (i) implies (ii) and that (iii) implies (i). 
We start with the latter. 

6.1. Proof that (iii) implies (i). By assumption, there is some p with £ — e < 9n < 6 
such that the stationary density has the form 

k 

p»(x)=J2 a i e ~ {Ci ' x) - ( 16 ) 
i=i 

We assume without loss of generality that the are nonzero and that the q are distinct. 
We may also restrict ourselves to the case of more than one summand (k > 1), since the 
k = 1 case has already been studied [T51 H7] . 

We next argue that p 11 satisfies (jl3|) — fjl5|) . To this end, in view of Proposition U we need 
to show that (BAR) is satisfied with di/Q = p^ do and dvi = p^/2dai. We do so using an 
argument due to Harrison and Williams [16, p. 108]. From Proposition [3] we know that 
(BAR) holds for some measures v\,V2- Let A £ R 2 satisfy (A,x) > for all x G S. On 
substituting f(x) = e~^ X ' x ^ in (BAR) we find 

i||A|| 2 + ( / u,A) [ e-^p»(x)dx-( Vl ,\) [ e- {x > x) vi(dx)-(v 2 ,\) ! e~^u 2 (dx)=0. 

2 J JS JF X J F 2 

Let us first focus on v\. Write A = Aii^o + A2^i, and let A2 — > 00 after dividing the above 
equation by A2. To evaluate the resulting limit of the first term, we use the initial value 
theorem to obtain 

lim A 2 / e- {x > x) p»{x)dx = / e - XlS p fl (sw Q )ds. 

After taking the limits of all other terms as well (recalling that (v\, n\) = 1), we conclude that 
J* Fi e~ XlXl v\{dx) = J °° e~ Xlz p^({z, Q))dz/2 for Ai > 0. The uniqueness theorem for Laplace 
transforms thus yields dv\ = p ti /2da\. A similar argument works to show dv 2 = p fi /2dd2 
by studying A = Xiw^ + A2«2 for large A2. 

In the remainder of this subsection, it is our aim to further specify the structure of p^ 
defined in (|16p when it solves (|13| ) — (|15[) . Lemmas [3HS] play a central role in this analysis. 

Graph representation. It is convenient to represent p^ by an undirected labelled graph 
G, with k vertices, as follows. Each vertex represents a summand aie~^ Ci ' x ^ in (|16p . and we 
say that two vertices are joined by a BC1 edge (respectively BC2 edge) if the sum of the 
terms corresponding to these vertices satisfies (fT3l) and (TH]) (respectively (fT3j) and (fT5l) ). 
Note that by Lemma O at most one BC1 edge and at most one BC2 edge can be incident 
at any given vertex. Therefore, the degree of the vertices in G cannot exceed two, and 
BC1 edges and BC2 edges alternate along any path. Here and throughout, we say that a 
subgraph of G is a path if it is connected and acyclic. The length of a path equals its number 
of vertices. The vertex corresponding to the summand e~^ c,x ^ is labelled by the matrix M 
for which (c, x) = (/i, (/ — M)x), and we refer to this vertex as an M vertex. Note that M 
is necessarily a reflection or rotation matrix. Although M is uniquely defined for any of the 
vertices of G, it has several representations — for example, we also refer to an vertex as a 
P2(Ro vertex. We refer to a path between an M\ vertex and and M2 vertex as an l M\ — M2 
path'. 

Mating procedure. Given the label of any vertex in G, we can specify the labels of all other 
vertices in the same connected component. For let Mi, M2 be the labels of two arbitrary 
vertices joined by a BC1 edge: then by LemmaO the sum of the corresponding exponential 
terms is proportional to p 7 for some 7 £ (0, 7r). We will say that Mi is the BC1 mate of M2. 
By considering separately the cases when M\ is a reflection and a rotation, it is easy to see 
that {Mi,M 2 } = {p2p, P2f3-2£,R{\ f° r some angle (3 £ (— tt,tt]. Similarly, the labels of two 
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BC2 BC1 BC2 BC1 

• • • m m 

P26 p2S^i P2S+2(, P25+2^R^ P2S+4£ 

Figure 3. An example labelled graph. By ([6]), this graph corresponds to 
the reflection construction in the leftmost diagram of Figure [2j If a = — 2 
then this graph also corresponds to the rightmost diagram of Figure [2 

arbitrary vertices joined by a BC2 edge are {p2/3, p2pR{\ f° r some j3 £ (— 7r, vr]. Any path in 
G — beginning for example with a BC2 edge — therefore has labels P2/3, p2pR%, P2/3+2£, • • • for 
some (3. 

Example. To illustrate the graph representation and mating procedure, suppose that one 
summand in (|16|) has exponent — (p, (I — p2s)x). In the graph representation, this summand 
is represented by a vertex with label p2S- suppose there exists a path of length 5 starting 
at this vertex, and that its first edge is a BC2 edge. By the mating procedure, the vertex 
labels for this path are p 2 s, P2sR^, P25+2§, P2S+2^R^, P25+4£, cf. Figure El By ©, this path 
corresponds to the reflection construction in the leftmost diagram of Figure [2l 

Proposition 5. Let the stationary density be of the form (|16|) . If G is the labelled graph 
corresponding to p^ then: 

• G is a P28 — p~2e path; 

• The number of vertices in G is odd; 

• a = — I for some integer i > 0. 

Proof. As usual we exclude the case a = 0. Since p^ is a density, each exponent in (]16p 
is nonpositive and p^ix) — > as — > 00 in 5. 

A key tool in the proof is the following range restriction for the reflection labels. That 
is, G cannot contain an vertex if w 7 € 5° or — w 7 £ 5°. Suppose a contrario that G 
contains an i? 7 vertex and that L 7 H S° 7^ 0, where L 7 is the line {rw-y : r G M}. We must 
then have p € L 7 since otherwise the exponent ih —(p,(I — R^)x) changes sign in S° on 
either side of L 7 . However, if p S L 7 then — (p, (I — i? 7 )x) = for all x, but such a constant 
exponent (which must be unique) contradicts p^{x) — > 0. 

G is acyclic. Suppose that G contains a cycle Go, of length 2m say (note that cycles 
of odd length are impossible by the mating procedure). Taking an arbitrary rotation label 
from Go, we must have p2p = P2/3+2m£ by cyclicity and the mating procedure, so that 
£ = mr/m for some integer n > 1. The reflection labels in Go are readily seen to be 
-Rg_l_fc£, k = 1, . . . , m, so range restriction yields Lp + ^ n S° = for each k = 1, . . . , m. Using 
this in conjunction with the fact that (by uniqueness of the reflection labels) {Lp + k n7T / m : 
k = 1, . . . , m} = {L j3+k7T / m : k = 1, . . . , m}, we deduce that n = 1 and that Rp+k^ = Ro for 
some k. By the proof of Lemma[5]i?o must be a vertex of degree 1, so this is a contradiction. 

G does not contain the label Ro, nor does it contain the label R^. We only prove that 
G does not contain the vertex Rq; similar arguments can be given for R^. Suppose that 
G contains the vertex Rq, and consider the sum-of-exponentials p^(x) corresponding to G 
when x lies on the boundary F\, There is a constant nonzero term since (/ — Rq)wq = 0, so 
to ensure p^(rwo) — > as r — > 00 there must be another exponent that vanishes on F\. First 
observe that Rq is the only reflection label that can possibly correspond to a constant term 
on Fi (as shown already, p L 7 if G contains the label R^ so p cannot be orthogonal to 
(/ — Ry)wo). All rotation vertices in G correspond to exactly the same exponents on F\ as 
their possible BC1 mates, which are reflection vertices necessarily different from Rq, whence 
none of the vertices joined by a BC1 edge can represent the term constant on Fx. Therefore, 
in view of Lemma U the only remaining possibility is that G contain a vertex labelled P2S 
which corresponds to the constant on Fx. However, this contradicts our assumption that 
£ - e < 0„ < S. 
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We have now proved that G is a p2S — p~2e path, and by the mating procedure the 
number of terms is odd, say 2t + 1. The mating procedure also shows that we must have 
P-2e = P2£(+28, so that 5 + e + ££, = wk for some integer n > 1. The range restriction on the 
reflection labels entails that L$+£, L$ +2 £, . . . , Lg + £g 5°. None of these lines can be equal 
to Lq, since this would imply that Ro is a label in G. We deduce that 5 + £(, < vr, and 
therefore n = 1. □ 

To continue our example, suppose that a = —2. Then /?2<5+4£ = P-2e and so the graph in 
our example also corresponds to the rightmost diagram in Figure EJ Then G has 5 vertices, 
so that has k = 5 exponential terms. 

6.2. Proof that (i) implies (ii). Suppose that a = — t for some t € {1,2,...}, and 
consider a p with 0^ G Qg. We shall use the representation (jHJ) and Lemma [H to argue that 
71"^ must equal the stationary density up to a multiplicative constant. We first argue that 
vr^ satisfies (TT3])~(fT5|) . 

The proofs of Lemmas 0] and [5] show that ir^ satisfies (I13D and (|14jl . To see that 7r^ also 
satisfies (j!5p we note that, for any constants do,. . . ,dg, the function 77^ defined on S by 

7r"(a0 = d (p, (I - p^e)v2)e-^ I - f, -^ x) + d 1 p i (x) + d 2 p 2§ (z) + . . . + d^(z) (17) 

satisfies (fl"3l) and ([Toll (pg is defined in Lemma EJ). Using e + <5 + ££ = 7r, on investigating the 
exponents we find that 7r^ and 7r^ are linear combinations of the same exponential terms, 
so it suffices to show that the coefficients are proportional to each other. To do so, we write 

tt"(x) = Q(M/-Rot^i)e-^ (/ - Rot ^ 

t 

+ [ci-kil*, (I ~ Rot^Me-^ 1 -^-^ 
fe=i 

- c M {p, (I - Ref^+O^e-^ 1 -^-^)*)" . 
Equating the coefficients with (|17j) . we find that 7r^ satisfies (|13p - (|15|) if for k = 1,...,£, 
C£^ k (p, (I-Rott-k)vi){pL, (I-Refe-k+ifa) = ci- k +i(p, (I-Reie-k+i)vi){p, (I-Rote- k )v 2 ), 

and this is the recursion given in Lemma [1] (the condition 9^ £ 0^ guarantees that none of 
the four inner products in ([10}) is zero). 

Now that we know that 71"^ satisfies (|13p ~(|15 p , it remains to show that it is a multiple of the 
stationary density. Proposition [2] and its proof show that 7r M is integrable and single signed. 
Moreover, Proposition Q] implies that 7r^(0) = 0. We therefore conclude from Proposition [J] 
that (BAR) is satisfied with dv$ = da and dvi = doi. The claim thus follows from 
Proposition [3j 

Acknowledgements 

We are grateful to Neil O'Connell for drawing our attention to (pQ), which motivated our 
research, and for several discussions. We would also like to thank Ivo Adan, Gerry Foschini, 
Mike Harrison, and Ruth Williams for stimulating discussions. This work was financially 
supported in part by the Science Foundation Ireland, grant number SFI04/RP1/I512. 

References 

[1] I. J.-B. F. Adan, J. Wessels, and W. H. M. Zijm, Analysis of the asymmetric shortest queue problem, 

Queueing Syst. 8 (1991), 1-58. 
[2] , A compensation approach for two-dimensional Markov processes, Adv. in Appl. Probab. 25 

(1993), 783-817. 

[3] F. Avrarn, J. G. Dai, and J. J. Hasenbein, Explicit solutions for variational problems in the quadrant, 

Queueing Syst. 37 (2001), 259-289. 
[4] P. Biane, P. Bougerol, and N. O'Connell, Littelmann paths and Brownian paths, Duke Math. J. 130 

(2005), 127-167. 



REFLECTED BROWNIAN MOTION IN A WEDGE 



13 



[5] J. Dai, Steady-state analysis of reflected Brownian motions: characterization, numerical methods and 
queueing applications, Ph.D. thesis, Stanford University, 1990. 

[6] J. G. Dai and J. M. Harrison, Reflected Brownian motion in an orthant: numerical methods for steady- 
state analysis, Ann. Appl. Probab. 2 (1992), 65-86. 

[7] A. B. Dieker and J. Warren, Determinantal transition kernels for some interacting particles on the line, 
Ann. Inst. H. Poincare Probab. Statist. 44 (2008), 1162-1172. 

[8] J. Dubedat, Reflected planar Brownian motions, intertwining relations and crossing probabilities, Ann. 
Inst. H. Poincare Probab. Statist. 40 (2004), 539-552. 

[9] P. Dupuis and K. Ramanan, A time-reversed representation for the tail probabilities of stationary re- 
flected Brownian motion, Stochastic Process. Appl. 98 (2002), 253-287. 
[10] P. Dupuis and R. J. Williams, Lyapunov functions for semimartingale reflecting Brownian motions, 

Ann. Probab. 22 (1994), 680-702. 
[11] G. J. Foschini, Equilibria for diffusion models of pairs of communicating computers — symmetric case, 

IEEE Trans. Inform. Theory 28 (1982), 273-284. 
[12] F. G. Foster, On the stochastic matrices associated with certain queuing processes, Ann. Math. Statistics 
24 (1953), 355-360. 

[13] J. M. Harrison, The diffusion approximation for tandem queues in heavy traffic, Adv. in Appl. Probab. 
10 (1978), 886-905. 

[14] J. M. Harrison and J. J. Hasenbein, Reflected Brownian motion in the quadrant: tail behavior of the 

stationary distribution, preprint, 2008. 
[15] J. M. Harrison and M. I. Reiman, On the distribution of multidimensional reflected Brownian motion, 

SIAM J. Appl. Math. 41 (1981), 345-361. 
[16] J. M. Harrison and R. J. Williams, Brownian models of open queueing networks with homogeneous 

customer populations, Stochastics 22 (1987), 77-115. 
[17] , Multidimensional reflected Brownian motions having exponential stationary distributions, Ann. 

Probab. 15 (1987), 115-137. 
[18] D. G. Hobson and L. C. G. Rogers, Recurrence and transience of reflecting Brownian motion in the 

quadrant, Math. Proc. Cambridge Philos. Soc. 113 (1993), 387-399. 
[19] W. Kager, Reflected Brownian motion in generic triangles and wedges, Stochastic Process. Appl. 117 

(2007), 539-549. 

[20] J.-F. Le Gall, Mouvement brownien, cones et processus stables, Probab. Theory Related Fields 76 (1987), 
587-627. 

[21] M. H. Protter and H. F. Weinberger, Maximum principles in differential equations, Prentice-Hall Inc., 

Englewood Cliffs, N.J., 1967. 
[22] S. R. S. Varadhan and R. J. Williams, Brownian motion in a wedge with oblique reflection, Comm. Pure 

Appl. Math. 38 (1985), 405-443. 
[23] R. J. Williams, Recurrence classification and invariant measure for reflected Brownian motion in a 

wedge, Ann. Probab. 13 (1985), 758-778. 
[24] , Reflected Brownian motion in a wedge: semimartingale property, Z. Wahrsch. Verw. Gebiete 

69 (1985), 161-176. 

[25] , Reflected Brownian motion with skew symmetric data in a polyhedral domain, Probab. Theory 

Related Fields 75 (1987), 459-485. 

University College Cork, Probability Group, 17 South Bank, Crosses Green, Cork, Ireland 

IBM Watson Research Center, Yorktown Heights, NY 10598 
E-mail address: ton.dieker@isye.gatech.edu 

University College Cork, Probability Group, 17 South Bank, Crosses Green, Cork, Ireland 



School of Mathematics, University of Manchester, Oxford Road, Manchester M13 9PL, 
UK 

E-mail address: j.moriarty@manchester.ac.uk 



